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Abstract 

Holographic interferometry and computer- 
assisted tomography (CAT) are used to determine 
the transonic velocity field of a model rotor 
blade in hover. A pulsed ruby laser recorded 
40 interferograms with a 2-ft-diam view field near 
the model rotor-blade tip operating at a tip Mach 
number of 0.90. After digitizing the interfero- 
grams and extracting fringe-order functions, the 
data are transferred to a CAT code.- The CAT code 
then calculates the perturbation velocity in sev- 
eral planes above the blade surface. The values 
from the holography-CAT method compare favorably 
with previously obtained numerical computations in 
most locations near the blade tip. The results 
demonstrate the technique's potential for three- 
dimensional transonic rotor flow studies. 


Notation 

a Q - speed of sound, ft/sec 

A - wave amplitude 

C » blade chord, ft 

f(x') - filter function 

I - irradiance 

k * Gladstone-Dale constant, ftVslug 

L * path length, ft 

n - refractive index 

n Q = ambient refractive index 
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N * fringe-order number 

R » spanwise coordinate, ft 

Rq - blade span, ft 

t - hologram amplitude transmittance 

U c - reconstruction wave complex amplitude 

Uj - transmitted wave complex amplitude 
U p - reference wave complex amplitude 
U Q1 » ambient object wave complex amplitude 
U 0 2 = test object wave complex amplitude 

Y - perturbation velocity, ft/sec 

x' - projection coordinate, ft 

X - chordwise coordinate, ft 

Y - height above blade centerline, ft 

S - film proportionality constant 

Y - ratio of specific heats 

A<J> - optical path-length difference (OPD) 

Q - field projection angle, deg 

X - laser wavelength, ft 

p 

p » air density, slug/ft 

p Q - ambient air density, slugs/ft^ 

n - blade rotational speed, rpm 



Introduction 

On many helicopters, the rotor blade’s 
advancing tip encounters transonic flow during 
forward flight. At these high Mach numbers, the 
rotor blade's performance suffers from compressi- 
bility effects that often cause shock waves to 
form near the blade tip; the shocks can extend to 
the acoustic far-field. Through theoretical and 
computational investigations, researchers attempt 
to understand the local shock generation of high- 
tip-speed rotors and its propagation to the far- 
field. However, because of the problem's complex- 
ity and the difficulty of obtaining detailed 
experimental information about the flow, accurate 
means for confirming transonic rotor-blade designs 
have been notably lacking. 

Shock waves cause a number of aerodynamic, 
dynamic, and acoustic problems on high-speed helir 
copter rotor blades. First, the shock rapidly 
increases the aerodynamic drag through energy 
dissipation, flow separation, and wave effects. 
Second, local shocks cause sudden large changes in 
pitching moment which can excite various blade 
torsional modes. As the blade rotates in forward 
flight, its Mach number and angle of, attack vary. 
The shock appears on the advancing side of the 
rotor disk and often results in large chordwise 
movements; these movements can be in opposite 
directions on the upper and lower surfaces as the 
Mach number and angle of attack change during each 
revolution. The changing shock positions on the 
upper and lower surfaces cause an unsteady load- 
ing, which produces fluctuating pitching moments. 
These moments can cause unexpected blade motions, 
oscillating loads on pitch links, and vibrations 
throughout the entire aircraft. Third, shock 
waves on an advancing-blade surface can "delocal- 
ize" 1 and extend directly to the far-field. Large 
amounts of acoustic energy radiate in front of the 
helicopter near the tip-path plane. This heli- 
copter impulsive noise is annoying in general and 
too easily detected in military applications. 

In an attempt to describe the transonic rotor 
flow field and to resolve the problems associated 
with it, promising theoretical models 2 and numeri- 
cal codes^' 1 * have been developed. The numerical 
codes compare favorably with blade-pressure mea- 
surements, 5 but are not yet verified at points 
away from the blade's surface. 

Previous attempts to measure the flow field 
have been severely limited. Pressure-instrumented 
airfoils are expensive and difficult to fabri- 
cate — especially in scale models. In addition, 
pressures can be measured only at the blade sur- 
face. Hot-wite anemometry requires that a probe 
be positioned within the field, therefore disturb- 
ing the flow. Laser velocimetry requires 


flow-seeding and, when shock waves are present, it 
is uncertain whether the seeds follow the flow 
faithfully. Both hot-wire anemometry and laser 
velocimetry can take only point-by-point measure- 
ments, requiring large amounts of running time to 
survey the rotor's three-dimensional field, a 
distinct disadvantage for rotor testing. 

Schlieren and shadowgraph photography provide only 
a qualitative two-dimensional representat ion of a 
three-dimensional flow. And Mach-Zehnder inter- 
ferometry provides quantitative information, but 
is extremely difficult to use in a large-scale 
experiment. Clearly, another experimental tech- 
nique which overcomes these limitations must be 
employed. 

Holographic interferometry is an effective 
diagnostic technique for making transonic^flow 
measurements . ^ Previous investigations^' in 
which two-dimensional flows over airfoils were 
studied show that accurate quantitative informa- 
tion is obtainable using holographic interferom- 
etry. However, the transonic flow around a heli- 
copter rotor blade is three-dimensional and 
requires a tomographic technique to compute the 
correct flow information from several interfero- 
grams . To date, most applications of this tech- 
nique have been limited to axisymmetric flow or to 
simple three-dimensional flow with a small model 

Q IQ 

under ideal laboratory conditions, ' 

This paper discusses the procedures necessary 
to obtain quantitative measurements of a tran- 
sonic, three-dimensional flow field near a rotor- 
blade tip, using holographic interferometric data 
and CAT. Though most helicopter rotor problems 
caused by shock waves occur during forward flight, 
this experiment investigates the steady problem 
(hover), which simulates many physical phenomena 
of forward flight. 11 The method for recording 
interferograms and example interferograms is 
included, and the steps that must be followed in 
extracting quantitative information from the 
interf erograms are outlined. The technique's 
potential for measuring three-dimensional tran- 
sonic rotor flows is demonstrated, and the results 
it yields are compared with those from previously 
performed numerical computations. 

Background Concepts 

For the experiment to be successful, it is 
necessary to 1) record high-quality interferograms 
near a rotor-blade tip from multiple viewing 
angles, and 2) implement a suitable CAT code with 
the interferogram data. Familiarity with holog- 
raphy, holographic interf erometry , and computer- 
assisted tomography principles provides the neces- 
sary insight for understanding this technique. 


2 



Holography 

Holography is a two-step imaging process in 
which diffracted light waves are recorded and 
reconstructed. 12,1 ^ The first step is recording, 
or storing, the hologram. This is accomplished by 
dividing a single coherent laser beam into two 
beams and exposing a photographic film to the two 
light waves, as shown in Fig. la. The object 
wave, which is the wave containing the flow infor- 
mation, passes through the, measured field (the air 
near the blade tip in this experiment). The sec- 
ond wave, the reference wave, passes around the 
field. By adding the coherent reference wave to 
the object wave, the photographic film records a 
high-frequency interference pattern. Once the 
film is developed, it is known as a hologram, 
which is a complicated diffraction grating. 

The second step in holography is reconstruc- . 
tion, or playing back the hologram. This is 
accomplished in two ways. First, a reconstruction 
wave identical to the reference wave illuminates 
the hologram (Fig. 1b). The hologram diffracts 
the reconstruction wave and produces a replica of 
the original object wave, forming the original 
object's virtual image. In the second method of 
reconstruction, a reconstruction wave conjugate to 
the reference wave illuminates the hologram 
(Fig. 1c). The hologram diffracts the conjugate 
wave forming the original object’s real image. 

The real image may be photographed without the use 
of a lens by placing a sheet of photographic film 
in the real image space. 

Several important characteristics of hologra- 
phy are applicable to the experiment at hand. 

There are very few geometrical constraints in a 
holographic optical system; thus, holography can 
be applied in a large-scale, nonlaboratory envi- 
ronment. Note that recording and reconstruction 
of the hologram can be done in different locations 
if the reference wave is reproducible . Thi3 
allows the reconstruction to be done in a labora- 
tory, far from the harsh environment in which the 
hologram was previously recorded. The reference 
wave serves only as a method of recording and 
reconstruct ing the object wave. Thus, a hologram 
does not produce quantitative information about 
the field of interest. To obtain quantitative 
information (in the form of interference fringes) 
an interf erogram must be recorded. 

Holographic Interferometry 

Holographic interf erometry is the interfero- 
metric comparison of two object waves recorded 
holographically (3ee the Appendix for further 
detail). In this experiment, the two object waves 
are recorded sequentially in time with double- 
exposure holographic interferometry. The inter- 


ferogram is recorded by first exposing a photo- 
graphic film to a reference wave and an "undis- 
turbed" object wave. Later in time, the same 
photographic film is exposed to a reference wave 
and to a second "disturbed" object wave. 

When the holographic interferogram i3 recon- 
structed, the virtual or real image shows the 
object (the transparent field) with an interfer- 
ence fringe pattern. The fringe pattern repre- 
sents the difference between the "undisturbed" and 
"disturbed" flow states. The irradiance of the 
reconstructed wave is proportional to 

1 - l U 01 + u o 2 | 2 

IP 10 

which can be written as * J 

I - 2A 2 [1 ♦ cos (Alt’) ] (1) 

Equation (1) represents the interferogram with a 
fringe pattern of dark and bright fringes of con- 
stant optical path-length difference (OPD) A<|>, 
where A<J> is given by 

A4> - J [n(x,y ,z) - n Q ] - NA (2) 

To determine the flow-field properties, the line 
integral of Eq. (2) must be inverted and solved 
for n(x,y,z), the refractive index at a specific 
point in the field. 

In a two-dimensional flow (i.e., the flow 
over a fixed airfoil in a wind tunnel), the evalu- 
ation of Eq. (2) is simplified, since the refrac- 
tive index is constant across the width of the 
test section L, Eq. (2) becomes 

n(x,y) - n Q ♦ NA/L (3) 

In a two-dimensional flow, the fringes on an 
interferogram are contours of constant refractive 
index and the refractive index at any point in the 
field can be determined from a single interfero- 
gram. However, since the transonic flow near a 
rotor-blade tip is three-dimensional, Eq. (3) 
cannot be used. To invert Eq. (2) and solve for 
the refractive index at a specific point in the 
field, computer-assisted tomography (CAT) must be 
employed. 

Computer-Assisted Tomography 

Tomography is a mathematical technique for 
reconstructing a three-dimensional field from its 
two-dimensional projections (see References 14 
and 15 for a wide variety of applications). A 
projection of a three-dimensional field is the 
fringe pattern recorded on an interferogram. All 
methods require multidirectional projection data 
of the field. Tomographic codes develop in two 
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directions: 1 ) iterative algebraic reconstruction 

techniques, 1 ® and 2) Fourier transform techniques. 
A version of the latter method, termed filtered 
back-projection, 1 ^' 1 ® appears most suitable for 
this application. 

Most Fourier transform techniques employ 
back-projection. Projection data from the field 
are recorded in one plane at several azimuthal 
angles around the field. For example, one projec- 
tion of a uniform absorbing disk is shown in 
Fig. 2a (taken from Reference 19). Beyond the 
disk boundary (no path length through the disk), 
the light ray’s OPD is unchanged, producing no 
interference fringes. Near the disk boundary 
(short path length through the disk), the light 
ray's OPD is changed slightly, producing a few 
interference fringes. And near the disk center 
(long path length through the absorbing disk), the 
light ray’s OPD is changed substantially, produc- ' 
ing several interference fringes. Similar projec- 
tions (fringe number vs position) at different 
azimuthal angles are also recorded. Each projec- 
tion is then back-projected, or smeared back along 
the direction in which it wa 3 recorded (Fig. 2b). 
Values are added point by point to form a recon- 
struction of the field. Unfortunately , simple 
back-projection produces an undesirable spoke 
pattern which severely degrades the quality of the 
reconstructed field. 

To eliminate the spoke pattern, the pro- 
jections are filtered. A one-dimensional con- 
volution (indicated by an asterisk) is performed 
between each projection and an appropriate filter 
function (see Reference 20 for a discussion of 
filter functions) before back-projection, as shown 
in Fig. 3a (taken from Reference 19). Each fil- 
tered projection is then back-projected over the 
reconstruction space (see Fig. 3b). The negative 
side-lobe 3 introduced by the filter eliminate the 
spoke pattern during the point-by-point addition 
process. With many projections, this technique 
yhields an accurate reconstruction of the original 
field. 


Procedure 

Several steps must be performed to quantita- 
tively reconstruct the three-dimensional transonic 
field near a model helicopter blade tip. First, 
several holographic interferograms must be 
recorded along planes perpendicular to the rotor 
tip-path plane at various azimuthal angles 
(Fig. *0. Data must then be extracted from the 
inter f erograms . This can be done 1) manually, 

2 ) by using a graphic tracing tablet , 21 or 3 ) py 
using a system that digitizes the interferograms 
and extracts fringe-order numbers. The digital 
interferogram evaluation technique was used; it 


will be presented in detail by Becker and Yu . 22 
Finally, the data are transferred as input to a 
tomography code, which computes the refractive 
index at specific points in a horizontal plane 
above the blade surface. This procedure is 
repeated in several planes to yield a reconstruc- 
tion of the entire three-dimensional field. 

Recording Holographic Interferograms 

The holographic system for recording inter- 
ferograms near a model rotor blade was assembled 
in the Aeromechanics Laboratory Anechoic Hover 
Chamber. Figure 5 shows a schematic of the opti- 
cal system and Fig. 6 shows the Anechoic Hover 
Chamber. A ruby laser with a 20-nsec pulse width, 
a 694.3-nm wavelength, and a power of 1 J 
’•freezes” the rotating blade at any desired azi- 
muthal angle. A beam-splitter divides the laser 
beam into two separate beams at the laser outlet. 

A microscope objective lens expands the object 
beam to fill a 24-in. -diam spherical mirror. 

Since the foci of both the objective lens and the 
spherical mirror coincide, a collimated plane wave 
forms as the beam passes through the rotor area. 
The object beam then strikes a second 24-in. -diam 
spherical mirror, emerges as a converging wave, 
and illuminates a 4-in. by 5-in. photographic 
plate. The reference beam is lengthened by caus- 
ing it to strike several plane mirrors. This beam 
must be lengthened so that the difference in the 
path lengths of the object and reference beams is 
les 3 than the coherence length of the laser (one 
of the very few, and easily met, geometrical con- 
straints in a holographic system). The reference 
beam is expanded by an objective lens, then colli- 
mated with a 5 -in . -dlam lens; finally, it is 
directed toward the film so that it overlaps the 
object beam. 

The entire procedure is conducted from out- 
side the hover chamber, once the optical system is 
aligned. Firing the laser, changing the photo- 
graphic film plates, and controlling the test 
conditions are all done by remote control. Recall 
that to record an interferogram , two exposures at 
different times (different flow states) must be 
made on a single film plate. The film records the 
first exposure while the rotor blade remains sta- 
tionary. In this case, the air has no velocity 
and therefore has a uniform refractive-index dis- 
tribution. The film records the second exposure 
while the blade rotates at the desired speed. The 
nonhomogeneous refractive-index distribution in 
this case introduces phase changes in the second 
object wave, producing interference fringes on the 
film plate. This double-exposure recording proce- 
dure repeats at various angles around the flow by 
synchronizing the laser pulse with the desired 
blade position. Because of the long optical path- 
lengths (90 ft), the recording system is very 



sensitive to vibrations of the optical components. 
At several azimuthal angles, it was necessary to 
record multiple interferograms to obtain one high- 
quality interferogram. The photographic plates 
are then removed from the recording system in the 
hover chamber, developed, and reconstructed in a 
laboratory for further processing. 

Holographic interferograms record the flow 
near a hovering 1/7-scale (geometric) model UH-1H 
rotor with untwisted NACA.0012 airfoil sections. 
The blade runs at a tip Mach number of 0.90 so 
that the flow is transonic and a shock wave is 
present. 1 This model normally uses two blades; 
however, in views along the span, the optical beam 
would pass through the refractive-index field of 
both blades. Because the refractive-index fields 
of the two blades are inseparable at these angles, 
a single-bladed rotor with a counterbalance is 
used instead (Fig. 7). 

Holographic interferograms near a transonic 
rotor blade are recorded at MO different viewing 
angles. The blade rotates in a clockwise direc- 
tion and can be captured at any desired viewing 
angle with the pulsed laser. The tomography code 
requires flow data from certain viewing angles 
within a 180° range. Numerical simulation 
results 2 ^ using numerical computations of the 
flow-^ suggest recording interferograms from 
0 =* 8° to 0 “ iJ0° and from 0 * 140° to 
0 = 186° in 2° increments, as defined in Fig. 8. 
The missing views, 0 - 42° to 0 » 138°, were 
presumed to have very few interference fringes and 
were not utilized. 

Illustrated in Fig. 8 are examples of holo- 
graphic interferograms recorded near the model 
blade tip. The fringe pattern’s appearance 
depends on the viewing direction. Interf erograms 
recorded along the chord (near 0 = 90°) display 
very few interference fringes, since the optical 
rays pass through the field’s thinnest (weakest) 
region. No observable details are present in 
these views. However, in views along the span 
(near 0-0° or 0 = 180°), numerous fringes are 
visible, because the optical rays pass through the 
longest (strongest) region within the field. The 
leading-edge stagnation point, shock structures, 
boundary-layer separation, and wake system are 
clearly visible. In particular, a lambda shock 
(0 * 180°) and the radiated shock (0 = 186°) 
appear above the blade. Several interferograms 
are described in detail in Reference 2H. 

Data Extraction 

One important step between recording the 
interferograms and applying tomographic recon- 
struction techniques is the evaluation of 
the interferograms; that is, reading fringe 


positions and fringe numbers. Up to now, most 
interferograms from aerodynamic tests have been 
evaluated by either reading fringe numbers and 
their positions manually or by tracing the fringe 
lines by hand with the help of a tracking device 
(for Instance a graphic tablet). Manual evalua- 
tion is a very time-consuming and inaccurate pro- 
cedure, however, and it is evident that in the 
current application, where large numbers of inter- 
ferograms have to be evaluated at several hori- 
zontal planes, utilization of an automatic fringe- 
reading procedure would enhance the evaluation and 
would make the interferometric technique a much 
more powerful measurement tool (see Reference 22 
for a detailed discussion of this procedure). 

In order to reconstruct the three-dimensional 
flow field from these interferograms, using tomo- 
graphic reconstruction techniques, the fringe 
number-functions along cross sections of a plane 
parallel to the rotor disk are required. Data 
from all the views are needed to reconstruct the 
index of refraction in a particular plane above 
the blade. Several planes have to be recon- 
structed to get a three-dimensional flow-field 
representation. This procedure requires that each 
interferogram be represented in a form that allows 
the fringe-order function to be computed at any 
desired point. 

The evaluation of interf erograms using com- 
puter-aided methods can be subdivided into the 
following steps: 1) digitization of the inter- 

ferograms and image enhancement, 2) fringe seg- 
mentation and fringe-coordinate extraction, 

3) merging of fringe fields obtained from several 
magnified views, *0 fringe numbering and correc- 
tion of fringe disconnections, 5) coordinate 
transformations , 6) interpolation and extrapola- 
tion of fringe-number functions. 

Hardware Components . An image-processing 
system (De Anza IP-61100) connected to a VAX 11/780 
host computer provides the main hardware necessary 
for digitizing the interferograms and for doing 
some image-enhancement processing. 

The resolution of the system is 512 x 512 
pixels with an 8-bit intensity range. Currently 
it is equipped with two of possible four memory 
planes, as well as with a graphic and an alpha- 
numeric overlay. A frame-grabbing unit can digi- 
tize a frame of a video signal in real time 
(1/30 sec). A black and white video camera (MTI 
series 68) with a resolution (bandwidth) of 18 MHz 
is connected to this input channel. 

The system has an arithmetic-logic unit 
(ALU), with which real-time addition, subtraction, 
or comparison of one or more image planes may be 
made. The contents of each memory plane may be 
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routed through lookup tables before being input to 
the ALU, to the video-output processor, or to 
another plane. The actual contents of any image 
plane or of a combination of image planes is 
output via a video-output processor and can be 
shown on a color-display. Each channel has its 
own color-mapping tables. A joystick control 
device is used for interactive input. It controls 
two cursors, which may be used in a number of 
different operating modes. A color print system 
(Dunn Instruments model 6}1 ) serves as a hardcopy 
device for the color monitor. 

Digitization and Processing . During the 
recording of the holographic interferograms , two 
fiducial points were marked in the image plane 
with a known position relative to the rotor 
system. This allowed a subsequent coordinate 
transf ormation into a space-fixed coordinate 
system. Each interf erogram was then digitized, 
and several subsections were enlarged, depending 
on the fringe density in the interf erograms . A 
scale, aligned to the fiducial marks, was always 
digitized together with the interf erogram so that 
the position and magnification of the enlarged 
segments could be identified. Before segmentation 
and f ringe-coordinate extraction are applied, it 
may be U3ef ul to do some image enhancement to 
reduce noise. Frame-averaging methods increase 
the 3 ignal-to-noise ratio of the imaging elec- 
tronic components, thereby improving the picture 
quality. 

One-Dimensional Evaluation . A one- 
dimensional interactive processing may be 
employed, if only one two-dimensional plane of the 
flow field is reconstructed; hence, only one cross 
section through each interf erogram has to be 
known. In this case some of the expense required 
for the full two-dimensional evaluation of the 
fringe pattern may be avoided. But the procedure 
has to be interactive, because local information 
(i.e., knowledge of the fringe locations along one 
line) is not sufficient to number the fringes 
correctly or to detect lost fringes or other dis- 
tortions. In an interactive procedure, informa- 
tion such a3 acceptance of the segmentation and 
assignment of fringe numbers has to be given by 
the user. 

A program featuring the one-dimensional eval- 
uation digitizes and preprocesses an interf erogram 
and does a fringe-segmentation along a line or 
along a set of lines through the field of view. 

The result of the segmentation procedure (a binary 
fringe pattern) is written back to the image 
screen for monitoring reasons (see Figs. 9a 
and 9b). The user may interactively change the 
acceptance threshold or edit the segmented cross 
section in heavily distorted regions before he 
continues to the numbering section. The cursor 


can be moved to each part of the segmented line, 
and fringe numbers may be assigned interactively 
by various commands. The task of assigning fringe 
numbers is supported by color coding the black 
parts of the fringes; this shows whether the 
fringe-order function increases or decreases by 1 , 
or if a discontinuity is present between adjacent 
fringes. The resulting fringe-order function may 
also be plotted onto a graphic terminal. A post- 
processing program merges all the data taken from 
different enlarged portions of an interf erogram 
(as shown in Fig3. 9a and 9b) and uses a spline- 
approximation to interpolate between the fringes. 
The final result is extracting the fringe order 
function at a desired height above the blade sur- 
face as shown in Fig. 9c. 

Two-Dimensional Interactive Processing . 

Often the knowledge of only a cross section across 
the field of view is not sufficient for interpret- 
ing a given flow problem; instead, the fringe 
patterns have to be evaluated over the whole field 
of view. To facilitate the two-dimensional inter- 
f erogram evaluation, additional algorithms have 
been implemented to segment fringes, to extract 
each fringe’s coordinates (trace the fringe 
sides), and to represent the whole fringe field as 
a polygonal data structure. Methods for numbering 
these line fields and detecting extraction errors 
have been developed. With a two-dimensional 
interpolation scheme, fractional fringe-order 
numbers may then be estimated at any given point 
in the field of view. 

Figure 10 shows an example of a typical 
fringe-coordinate extraction. An interf erogram 
(Fig. 10a) is digitized over the whole field of 
view with low resolution. The white lines over- 
laying the display represent the extracted fringe 
sides, and the black line represents the border 
polygon. The border polygon was input by hand 
before the coordinate extraction in order to mask 
off background objects and the part of the fringe 
field with high fringe density. (In our implemen- 
tation, all the graphics are In color and there- 
fore easier to distinguish from the displayed 
Image. ) 

As was discussed before, it may be necessary 
to digitize enlarged subsections to resolve 
regions with high fringe densities. The resulting 
fringe polygonal field is then combined out of all 
the individual digitized and processed subsec- 
tions, starting with the one with highest resolu- 
tion. To handle the problem of partly overlapping 
polygonal line fields, a boundary polygon is main- 
tained for each subsection. This defines the 
definition area for each fringe polygonal field. 
Upon combination of two overlapping fields, the 
border polygon of the one with higher resolution 
(priority) is used to intersect the lines of the 
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other field. Both line fields are then connected 
at those intersection points. The new border 
polygon is the border of the combined areas. 
Figure 10b shows an enlarged digitized part of 
the interferogram shown in Fig. 10a, with the 
extracted fringe polygons overlaying. 

CAT Reconstruction 


Fringe-order functions are transferred to the 
filtered back-projection CAT code which computes 
the refractive-index field at specific points in a 
chosen horizontal plane above the blade surface. 
The code assumes refractionless light rays; there- 
fore, each horizontal plane can be treated inde- 
pendently, even though data for each plane is 
taken from one set of interferograms. The pertur- 
bation velocity is computed by first converting 
refractive index to density, using the Gladston-^ 
Dale relation: 


p - (n - 1 )/k 


Density is then converted to perturbation velocity 
from a form of Bernoulli's equation for steady 
(with respect to the rotation blade), compressi- 
ble, isentropic flow: 



The procedure is repeated in several planes above 
the blade to reconstruct the entire three- 
dimensional field near the model blade tip. 


Reconstruction Results 


The holography-CAT reconstruction of the 
blade-tip velocity field is compared with numeri- 
cal computations. The computations used here are 
conservative, mixed-difference solutions of the 
transonic small-disturbance equation. Results 
from both sources are presented in four horizontal 
planes, as identified in Fig. 11a. Three plot 
types are used to visualize the flow field. 

First, velocity contours are given in plan view 
(see Fig. 11b), where the blade's leading and 
trailing edges are at X/C « 0.0 and X/C - 1.0, 
respectively. The blade tip is located at 
R/Rq « 1.0, the rotation center is at R/R 0 " 
and the blade rotates in a clockwise direction. 
Second, perspective views are displayed in which 
velocity values are plotted along the vertical 
axi 3 . The data and geometry are identical to the 
contour plots, though the data are viewed from 
near the rotor hub. Third, velocity distributions 
are shown at six radial locations (see Fig. 11b) 
for each plane. 


Figure 12 compares the velocity contours 
derived from the holographic-CAT method and the 
numerical computations near the blade surface 
(Y/C - 0.08). Both methods display low-velocity 
regions near the leading and trailing edges, and 
also display a high-velocity region over the blade 
surface containing a shock at approximately 
X/C * 0.60 near the blade tip. The general con- 
tour shapes show a strong resemblance except near 
the blade tip (roughly the last 556 of blade span) . 
The maximum velocity region appears at the shock 
foot in the computational analysis, but it appears 
closer to the leading edge and farther from the 
shock foot in the holography-CAT results. 

Figure 13 shows velocity values for the same 
plane in perspective view. Again, the general 
flow shapes appear very similar. The major dif- 
ference between the two results is the roughness 
(minor "ridges") in the reconstructed flow. This 
may be due to reconstruction artifacts caused by 
noncontinuous data (interferograms recorded in 2° 
azimuthal increments) or by noise (erroneous 
fringes) in the interferogram data caused by opti- 
cal component motion. Figure 14 compares velocity 
distributions at six radial locations. The rough- 
ness of the reconstructed flow can be seen 
throughout the figures. Also, the holography-CAT 
method determines the shock location to be 
slightly more downstream (350 than does the numer- 
ical code solution. The major differences can be 
observed in Fig. I4d, where the discrepancies at 
the leading edge and over the blade surface are 
clearly visible. The leading edge (X/C =* 0.0) 
difference may result from a breakdown of the 
small crossflow assumption [Eq. (4)] near the 
blade tip. The difference over the blade surface 
(X/C - 0.2 to X/C - 0.6) may be attributed to 
the existence of a lambda shock, in the interfero- 
gram data (i.e., Fig. 8f), which cannot be pre- 
dicted by the nonviscous numerical potential code. 
A lambda shock was also observed in Schlieren 
photographs from a previous wind-tunnel test 
(Fig. 15, taken from Reference 25) using the same 
airfoil and tip Mach number. 

Figure 16 compares the velocity contours from 
both the holography-CAT method and numerical code 
at Y/C - 0.22 (near the upper region of the 
lambda shock). The general velocity contours show 
an excellent agreement in both shape and magnitude 
throughout the plane. The maximum velocity region 
on the blade surface (near X/C - 0.50, 

R/Rq - 0.96) match much closer in this plane than 
in the plane near the blade surface (Fig. 12). 

The perspective view of Fig. 17 also shows an 
excellent agreement in shape and magnitude. 

Again, the most noticeable difference is the extra 
ridges in the reconstructed flow at the 3ame loca- 
tions and orientations as seen in the previous 
plane (Fig. 13a). The velocity distributions in 
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this plane (Fig. 18) compare favorably, especially 
inboard of the blade tip (see Figs. I8a-l8c) t as 
well as at the blade tip (Fig. I8d). The lambda 
shock's effect near the tip is apparently weaker 
in this plane; thus, the velocity distribution 
magnitudes and shapes are much closer than those 
in Fig. 14. 

Figure 19 compares velocity contours at 
Y/C = 0.49 above the blade. The velocity contour 
shapes and magnitudes are similar in all regions 
except that the holography-CAT method shows the 
maximum velocity point to be slightly span) 
outboard of the numerical result. The extra 
ridges seen in the perspective view (Fig. 20) are 
in approximately the same location and orientation 
as those in the previous plane; however, the mag- 
nitude of the ridges has decreased. Figure 21 
shows velocity distributions for this plane. 

There is good agreement at all locations except * 
beyond the blade tip (Fig. 21e) where the holog- 
raphy-CAT method shows slightly larger velocity 
magnitudes . 

Finally, Fig. 22 shows velocity contours 
for Y/C = 1.17. The velocity contours match 
throughout, though there are no distinguishing 
features in the flow at this height above the 
blade surface. In this plane, the extra ridges 
are almost unobservable in the perspective view 
(Fig. 23). The velocity distributions of Fig. 24 
also show a 3trong similarity between the two 
methods at all radial locations, for both results 
show that this plane is at the perturbed flow's 
upper extent. 

Overall, the agreement between the holog- 
raphy-CAT results and the numerical solution is 
extremely encouraging. However, comparisons with 
other experimental data sources are required 
before a final decision can be made about the 
holography-CAT results. Several discrepancies 
must be resolved. First, data from pressure- 
instrumented blades will aid in confirming the 
shock location and whether a lambda shock does 
exist in the flow. Second, the extra ridge3 
appearing in the reconstructed results must be 
eliminated or reduced; the ridges may be caused by 
poor data recorded in the interf erograms. To 
improve the quality of the i nterf erograms , a modi- 
fied optical system (for both hover and forward 
flight testing) is necessary. In addition, 
neither the holography-CAT reconstruction nor the 
numerical code solutions compute the expected 
shock strength beyond the blade tip. Acoustic 
measurements indicate a much stronger radiated 
shock than is indicated by these two results. 


Concluding Remarks 


The holographic i nterf erometry computer- 
assisted tomography technique proved to be a 
highly effective way of measuring the three- 
dimensional, transonic flow field near a model 
rotor-blade tip. Results from this method compare 
favorably with those of numerical computations, 
except very near the tip region. That discrepancy 
may be due to the existence of a lambda shock 
recorded by the interf erograms (and in previous 
Schlieren photographs) which is not predicted by 
the nonviscous potential code. In other regions, 
the velocity distributions along the chord are 
similar in both shape and magnitude. However, the 
results from the technique must be further' veri- 
fied against other experimental data. 

Since this is the first successful implemen- 
tation of the holographic interferometry , com- 
puter-assisted tomography method in rotor flow 
studies, many improvements are indicated. For 
example, the optical system must be improved so 
that better quality 1 nterf erograms can be 
recorded, and an automatic fringe-reading tech- 
nique must be completed so that the time required 
to evaluate interf erograms can be shortened. Upon 
verification of these results and after the system 
is improved, measurements of other model rotor- 
blade flow fields, including those of forward 
flight, will be performed. 


Appendix: Holographic Interf erometry 

Double-exposure holographic interf erometry is 
the interf erometric comparison of two object waves 
that are recorded holographically sequentially in 
time. The interferogram is recorded by first 
exposing a photographic film to a reference wave 
and an "undisturbed” object wave, as shown in 
Fig. 25a. Later in time, the same photographic 
plate is exposed to a reference beam and a second 
"disturbed” object wave, as shown in Fig. 25b. 

The irradiance at the film plane is 


K + Vl + K + VI 


2I r + V + V + W + V> 


+ VV + V* 


( A1 ) 


The last two terms of Eq. (A1) represent the . 
interference pattern recorded on the film, which 
contains both amplitude and phase information 
about the two reference and the two object waves. 
The amplitude transmittance of the developed film, 
called a holographic interf erogram, is 
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■ < 2I r * I 0, * V + U r ,(U 01 + U 02> 


+ V U 01 + V* 


When the reconstruction wave illuminates the 
interferogram, a3 shown in Fig. 25c, the trans- 
mitted light is 


U e fc 


’ U c (2I r + I 0, + V * U ? U c (U 01 + V 


U U (U„, 
r c 01 


U )* 
02 ; 


(A2) 


The second term in Eq. (A2) i3 the "composite” 
virtual image reconstruction of the two object 
waves (U Q1 ♦ U Q2 ). By illuminating the inter- 
ferogram with a conjugate reconstruction wave 
(Fig. 25d), the transmitted light is 


U*t 

c 


U*(2I + 

c r 01 


I 02 ) 


W U 0, + U 02 ) 


W U 01 


U )* 
02 ' 


(A3) 


The third term in Eq. (A3) is the composite real 
image reconstruction of the two object waves 
(Uqi + U Q2 ), which is the image that is photo- 
graphed and presented in this paper. 


The primary reason for using holographic 
interferometry Is that it possesses a property 
called cancellation of path-length errors. In a 
holographic interferometer , the interfering waves 
divide temporally. That is, the two interfering 
(object) waves are recorded at two different times 
but travel the same path through the optical sys- 
tem (recall that the reference wave serves only to 
3tore and play back the interf erogram) . Since 
there is no difference between the two path 
lengths of the interfering object waves in the 
optical system, only changes in the path lengths 
caused by the different states of the test field 
are displayed as fringes in the interferogram. 

This leaves a useful interferogram, even though 
low-quality optics are used. Therefore, holo- 
graphic interferometry can be used in a large- 
scale experiment in which relatively low-quality 
optics are used and still yield high-quality 
interferograms. 
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Fig. 1 Optical ho! 


Fig. 2 Back-pro jec 
each projection 






Fig. 5 Schematic drawing of the holographic recording system. 



Fig. 6 Holographic setup at Anechoic Hover 
Chamber. 


Fig. 7 One-blade rotor with a counterweight 
balance. 
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a) b) 


HEIGHT 0.673, ANGLE 182.50, M = 0.900 



Fig. 9 One-dimensional interf erogram evaluation, a) Digitized interf erogram with segmented cross 
section superimposed; b) evaluation of an enlarged section of the interf erogram shown in Fig. 9a; 
c) fringe-order function along cross section 3 after combination of several enlarged subsections and 
rational spline interpolation. 
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Fig. 10 Example of fringe polygon extraction, 
a) Original digitized i. nterf erogram — extracted 
fringe polygons are overlayed in white, the border 
polygon is written in black; b) enlarged digitized 
part of the same interf erogram — extracted fringe 
polygons are overlayed. 


Fig. 11 Rotor flow-field geometry, a) Location 
of four horizontal planes for velocity contour 
maps and perspective views; b) location of six 
radial stations for velocity distribution plots. 





Fig. 12 Perturbation velocity contours in plan view for Y/C - 0.08 above blade centerline: contour 

interval - 50 ft/sec. a) Holography-CAT reconstruction; b) numerical code. 



(a) 


(b) 


Fig. 13 Perturbation velocity values for Y/C - 0.08 above blade centerline, a) Holography-CAT 

reconstruction; b) numerical solution. 












Fig. 16 Perturbation velocity contours in plan view for Y/C = 0.22 above blade centerline: contour 

interval ~ 50 ft/sec. a) Holography- CAT reconstruction; b) numerical code. 
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Fig. 19 Perturbation velocity contours in plan view for Y/C = 0.49 above blade centerline: contour 

interval = 35 ft/sec. a) Holography-CAT reconstruction; b) numerical code. 




Fig. 20 Perturbation velocity values for Y/C - 0.49 above blade centerline, a) Holography-CAT 

reconstruction: b) numerical solution. 






x/c x/c 

Fig. 21 Perturbation velocity distributions at six radial locations for Y/C - 0.49 above blade 
centerline, a) R/Rq - 0.88; b) R/R Q - 0.92; c) R/R Q = 0.96; d) R/R 0 - 1.00; e) R/R 0 = 1.04; 
f) R/R 0 » 1 .08. 
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Fig. 23 Perturbation velocity values for Y/C * 1.17 above blade centerline, a) Holography-CAT 

reconstruction; b) numerical solution. 












Fig. 25 Double-exposure holographic interferometry, a) First-exposure recording; b) second-exposure 
recording; c) reconstruction of the virtual image; d) reconstruction of the real image. 
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